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Abstract. Dark matter haloes in cosmological N-body simulations are affected by processes 
such as mergers, accretion and the gravitational interaction with baryonic matter. Typically 
the analysis of dark matter haloes is performed in spherical or elliptical bins and the velocity 
distributions are often assumed to be constant within those bins. However, the velocity 
anisotropy, which describes differences between the radial and tangential velocity dispersion, 
has recently been show to have a strong dependence on direction in the triaxial halos formed 
in cosmological simulations. In this study we derive properties of particles in cones parallel or 
perpendicular to the collision axis of merger remnants. We find that the velocity anisotropy 
has a strong dependence on direction. The finding that the direction-dependence of the 
velocity anisotropy of a halo depends on the merger history, explains the existence of such 
trends in cosmological simulations. It also explains why a large diversity is seen in the velocity 
anisotropy profiles in the outer parts of high-resolution simulations of cosmological haloes. 
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1 Introduction 

In the ACDM-paradigm the first bound dark matter haloes are formed through the collapse 
of overdense regions, which have decoupled from the Hubble expansion. These haloes grow 
when they accrete matter from their surroundings and when they merge. Understanding 
merger remnants is therefore essential in order to understand the dark matter distribution 
in the universe. 

The dark matter (DM), making up most of the matter in the universe [1-3], consists of 
particles interacting so weakly, that haloes can be well modelled with collisionless mechanics. 
An example of evidence for this collisionless nature of DM comes from observations and 
simulations of the Bullet Cluster [4-7], which consists of two merging galaxy clusters with 
the collisional baryonic matter separated from the collisionless dark matter. 

A difference between collisional and collisionless particles is that collisionless parti- 
cles tend to have anisotropic velocity distributions, typically parametrised by the velocity 
anisotropy parameter [8], 

2 <d( r ) 
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where <7t a n is the total tangential velocity dispersion at a given radius (r), and cr ra d is the 
radial velocity dispersion. The velocity distributions are radially dominated for /3 > 0, 
tangentially dominated for (3 < 0, and isotropic for (3 = 0. In cosmological simulations the 
/3-profiles increase from f3 ~ in the central regions to (3 ~ 0.25 at r_2 [9, 10], which is the 
radius with a density slope of 7 = —2, where 

7=^. (1.2) 

alogr 

At larger radii the /3-profiles vary from halo to halo. 

Dark matter haloes from cosmological simulations have universal pseudo-phase-space 
density profiles following power laws [11], 

p/o^r) oc r~ a , (1.3) 

with a ~ 1.91 [10]. The density profiles are also universal [12, 13], with 7 ~ — 1 in the inner 
parts and 7 ~ —3 in the outer parts. Typically haloes are parametrised with the NFW-profile, 
or the Einasto profile [14, 15]. The universality of the density- and the pseudo-phase-space- 
density-profile in cosmological haloes determines cfad( r )- The velocity anisotropy, however, 
remains undetermined. 

In an attempt to understand the physics of the velocity anisotropy, several studies have 
examined /3-profiles in non-cosmological simulations. A linear relation [16, 17] between /3 and 
7 is consistent with the inner parts of cosmological haloes. More recently a 1-dimensional 
relation (an attractor) between /3, 7 and the slope of the radial velocity dispersion profile 
has been found in controlled non-cosmological simulations [18, 19]. A tight relation between 
/3 and 7, which is predicted by this attractor, is, however, not consistent with the large 
scatter seen in the outer parts of cosmological haloes. In the two proposed relations, which 
determines /3(r) from a given density profile, 7(7"), it is implicitly assumed that (3 is constant 
in spherical or elliptical bins. 

In a detailed study of a cosmological halo [20], it has, however, been found that /3 
behaves differently along different axes. The aim of this paper is to show that such direction- 
dependent /3-profiles can arise because of mergers. We will address this issue by studying the 
velocity anisotropy in different directions of merger remnants, to see whether /3-profiles are 
direction dependent. In our analysis we will cut out cones of particles in different directions, 
and then calculate (3 and other variables for the particles inside each cone. In Section 2 
we will present merger simulations, and in Section 3 we will map the velocity anisotropy 
parallel and perpendicular to the collision axis and determine whether they are different. 
We also show how structures are affected by smooth accretion (Section 4), and skymaps of 
the /3-profiles are presented (Section 5). Section 6 discusses relations between f3 and 7, and 
Section 7 examines pseudo-phase-space density profiles. We also discuss the shape of the 
merger remnants (Section 8). 

Throughout this paper we will use units with G = 1 and let log x denote the logarithm 
with base 10. 

2 The merger simulations 
2.1 Initial conditions 

We are interested in doing a simulation of the merging of two identical dark matter haloes. 
To do so Eddington's formula [21] (with /3 = 0) is used to generate a halo with a Hernquist 
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density profile [22], 



p (r) = -P* (2.1) 

PV ; r/r s (1 + r/r s ) 3 v 7 

We choose po = 1/2-k and r s = 1, which gives a structure with a total mass of 1. The density 
profile is truncated at a radius of 10r s . This density profile is chosen because of its simple 
analytical properties, and because it resembles [23] the NFW-profile, which describes haloes 
in cosmological simulations. 

As initial conditions two such structures are generated, and placed 20r s away from each 
other along the x-axis. 10 6 collisionless particles are used to represent each halo. The haloes 
approached each other with an initial velocity of 80% the escape velocity at this distance for 
an isolated Hernquist halo. Simulations were run with no impact parameter, and with an 
impact parameter of 10r s (in the y-direction) and with the same initial velocity and initial 
offset along the x-axis. 

Simulations were also made with an Osipkov-Merritt model [24, 25] (still following the 
Hernquist profile described above) with the velocity anisotropy given by 

2 

P(r) = (2.2) 
r A + r| 

where r a is the anisotropy radius, which we fix to 1.2r s . An anisotropy radius slightly larger 
than r s is chosen to avoid instabilities [26]. 

Simulations with 1:10 mergers were also performed (only for the model with (3 = 0). A 
structure, identical to the (3 = model defined above, was collided with a Hernquist structure 
with a 10 times lower total mass and a r s -value and cutoff radius 3 times smaller. 

2.2 Simulation details 

The public version of the N-body simulation code Gadget-2 [27] was used to run the simula- 
tions. All particles were collisionless and a spline softening of 0.015 was used in all simulations. 
All simulations were run for 300 time units. 

3 Analysing particles in cones 

To map the velocity anisotropy in the remnants, we defined cones pointing along the x-, y- 
and z-axis. The angle between the position vector (r) of a particle and a unit vector (n) 
pointing in the direction of a cone is 9 = arccos(n • r/|r|). We chose the cones to have an 
apex angle of 45°, so a particle is inside a cone if 9 < 22.5°. The cones are centered in the 
center of the merger remnants, i.e. at the position of the particle with the lowest potential. 
Figure 1 shows the particles in two cones for one of the merger remnants. 

3.1 Major mergers 

First the remnant of the major merger simulation for the /^initial = 0-haloes without an 
impact parameter is studied. Figure 2 (left panel) shows the velocity anisotropy in cones 
along each axis. Also shown is the spherically averaged value of (3(r). In the cone pointing 
in the direction of the collision axis (x-axis) the /3-profile is clearly different from (3 in the 
two perpendicular cones along the y- and z-axis. In the inner parts (with r < 1) (3 is almost 
constant ((3 ~ 0.2 — 0.4) along the collision axis, whereas it is an increasing function of radius 
in the two other cones. 
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Figure 1. Properties of particles inside a cone parallel and a cone perpendicular to the collision axis 
are studied. Particles outside the cones are not shown. 

In the figure an arrow marks the radius of r_2, which has been found by examining j(r). 
For a Hernquist profile, which was used as initial conditions in the simulation, r_2 = r B /2. 
In cosmological simulations the virial radius is typically 10r_2 [28], and in our simulations 
this is also roughly the radius, to which structures are equilibrated. 

The radial (cr ra d)> tangential (<7 tan ) and total (a) velocity dispersions for each of the 
three cones are shown in the right panel of Figure 2. Most striking is the significantly lower 
value of the tangential dispersion along the collision axis. The radial dispersion profile is 
nearly the same in all the cones. 

Figure 3 shows how the velocity anisotropy and the velocity dispersions (upper panels) 
depend on the local density in each cone. Using p as the abscissa corresponds to binning 
particles of equal density. f3(p) and o~ 2 (p) are clearly different in the three cones, so the 
/3-profiles and the velocity dispersions are not constant along contours of equal density. 

Another observation (lower panels) is that the slopes, 7(7"), of the density profiles are 
similar along the different axes, but the actual value of p, is significantly larger in direction 
of the collision axis. 

3.1.1 Comparison with the Via Lactea II halo 

In the Via Lactea II simulation [29] local properties of a Milky- Way-like halo have been 
examined [20]. A direction-dependent /3-profile is found in the inner parts (at 8 kpc) of the 
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Figure 2. Properties of particles in cones along the x-, y- and z- axis are analysed together with the 
particles in spherical bins. The merger has no impact parameter, and the collision axis is the x-axis. 
Anitiai = for the haloes. The left panel shows the velocity anisotropy in each cone, and the right 
panel shows the total, the tangential and the radial velocity dispersion in each cone. The grey arrow 
shows the radius, where 7 = — 2. 

halo: along the major axis f3 is larger (0.2 < (3 < 0.45), than along the intermediate axis 
(-0.1 < p < 0.25) and the minor axis (-0.65 <P< -0.1). 

Our merger remnant is elliptical with a major axis pointing along the ^-direction (see 
Section 8 for details about the shapes). We can therefore confirm the finding (from [20]) that 
(3 is largest along the major axis. 

3.1.2 The redistribution of kinetic energy 

Imagine two merging haloes, A and B, approaching each other with a velocity similar to 
the escape velocity at the cluster's separation distance. Seen from the center of mass of 
halo A, the particles from halo B have much larger kinetic energies along the collision axis 
than along the perpendicular axes due to the relative motion of the haloes. The reason why 
merger remnants have different /3-profiles in different cones is likely that merging processes 
are inefficient at transferring the kinetic energy along the collision axis, into kinetic energy 
along the other axes. Many particles in a merger remnant are likely left on non-spherical 
orbits, where they oscillate back and forth along the collision axis. Such a scenario is perfectly 
consistent with what we see in Figure 2. 

To test this explanation further the kinetic energies (K x , K y and K z ) along each axis 
are calculated in each cone, see Figure 4. Along the collision axis (x), we find that K y ~ K z 
at all radii. Along the y-axis, K x is larger than K z at all radii. This supports our claim, that 
the kinetic energy is not totally redistributed from the direction of the collision axis to the 
perpendicular directions in a merger. 

An observation is that K x ~ K y in the outer parts of the cone along the y-axis. This 
means that in a cone perpendicular to the collision axis, we have that <r(collision axis) ~ er ra d 
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Figure 3. /3(p), <J 2 (p), p(r) and j(r) in the different cones for a major merger without an impact 
parameter. 



in the outer parts. 

3.2 The effect of an impact parameter 

The presence of impact parameters breaks the rotation symmetry around the collision axis 
(which was present in the simulation without impact parameters), since the merging haloes 
are spiralling around each other in the x-y plane. 

Figure 5 shows the anisotropy and the velocity dispersions, when an impact parameter 
of 10r s is present. The rotation velocity at a given radius was subtracted before the velocity 
dispersion in a given bin was calculated. The velocity anisotropy is constant and positive in 
the inner regions in both the x- and in the y-direction, and negative in the z-direction. It 
is different from the simulations without impact parameters, where j3 was negative in both 
the y- and the z-direction. The spherically averaged /3-profile has a non-monotonic shape, 
similar to what is found in other studies of major merger remnants (e.g. [30]). 

The velocity dispersions are similar in all the three cones, except for the inner parts, 
where the cone in the z-direction has a larger tangential and a lower radial dispersion. 
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Figure 4. The mean kinetic energies in the x-, y- and z-direction along the three axes. K x in the 
cone along the y-axis is larger than K y in the cone along the x-axis. This shows that there is more 
kinetic energy along the collision axis of the merger remnant than along the other axes. 




Figure 5. The velocity anisotropy and the velocity dispersion for the major merger simulation with 
an impact parameter. The merging process is taking place in the x-y plane. 



3.3 Major mergers with Osipkov-Merritt haloes 

Figure 6 shows the remnant of two Osipkov-Merritt haloes collided without an impact param- 
eter. In this case the merger remnant again had different velocity anisotropies and velocity 
dispersions in the different cones. 

Before the progenitor cores touched each other, the two haloes were clearly more elon- 
gated, when visually inspected, than for the Anitial = mergers. Instabilities, such as the 
radial orbit instability, are therefore likely play an important role, and such instabilities might 
be the reason why the velocity anisotropies are not the same in the y- and the z-direction, 
which would be expected from the symmetries of the initial conditions. 
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Figure 6. Same as Figure 2, but with initial haloes following Osipkov-Merritt models. 



3.4 Minor mergers 

Figure 7 shows the velocity anisotropies in a 1:10 merger remnant for a simulations without an 
impact parameter. Three sets of particles are analysed: i) all the particles in the simulation 
(left panel), ii) the particles that started in the small halo (central panel) and iii) the particles 
originating from the main halo (right panel). The particles from the main halo remain at 
(3 = 0, and are almost unaffected by the merging process. The particles from the small halo 
are in orbits with f3 = 0.3 — 1.0. When all particles are analysed f3 is positive ((5 ~0.4) along 
the collision axis, but a nearly isotropic velocity distribution (f3 ~ 0.1) are present in the 
perpendicular directions. 

With an impact parameter present the velocity anisotropies are as shown in Figure 8. 
The main halo particles remain at ft = 0, but the spherically averaged /3-profiles are again 
positive due to the contribution from the particles from the small halo. Also note that the 
velocity anisotropy in the y-direction is larger than in the simulation without an impact 
parameter. 

The complicated direction-dependence of the velocity distributions in the minor merger 
remnants, is likely similar to the direction-dependent velocity structure in the substructure 
of cosmological simulations [31]. 

4 The effect of an accretion process 

Like mergers, smooth accretion of matter is a very important process for the build up of 
structure in the universe. We now want to study how the asymmetric /3-profiles in merger 
remnants are affected by a spherically symmetric perturbation that mimics smooth accretion. 

In the remnant of the major merger simulation without an impact parameter (from 
Section 3.1), we inserted particles in a spherical shell in the region, 4.8 < r/r_2 < 15, with 
a density profile following, p(r) oc r~ 2 . Each particle has a radial infall velocity identical 
to the escape velocity at its position. The number of particles in the shell and the mass of 
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Figure 7. Particles in a remnant of a 1:10 merger without an impact parameter are analysed in 
cones. The left panel shows all the particles in the simulation, the central panel shows particles that 
started in the small halo, and the right panel shows particles from the main halo. 
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Figure 8. Same as Figure 7, but with an impact parameter along the y-axis. 

each particle are the same as in the major merger simulations (2 x 10 6 particles with a total 
mass of 2). The simulations were run for 300 time units, which corresponds to 250 dynamical 
times, y/r 3 /GM(r), at r_2 of the final structure, and 20 dynamical times at 10r_2- 

Figure 9 shows the evolution of the density profiles for the particles that started in the 
merger remnant, and for the infall particles that started in the shell. The infall particles only 
gives a minor contribution to the central densities of the final structure (< 5%). 

Figure 10 shows the anisotropy for the final structure for all the particles, the accreted 
particles and the particles that started in the merger remnant. For all sets of particles, j3 is 
largest in the x-direction, and smaller in the directions perpendicular to the merger axis. An 
interesting observation is that the /3-profiles of the accreted particles are very similar to the 
/3-profiles of the particles from the merger remnant. 
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Figure 9. The initial conditions and the results of the simulations, where a merger remnant (grey 
points) is perturbed by a process that mimics smooth accretion. A group of particles accretes onto 
the structure from a shell at radii 4.8 < r/r_2 < 15 (black circles). 
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Figure 10. The velocity anisotropy of the particles in the accretion simulation. The left panel shows 
P(r) for all the particles in the simulation, the central panel shows the particles from the accretion 
shell, and the right panel shows the particles that started in the merger remnant. 



We conclude that the asymmetric velocity anisotropies remain asymmetric, when per- 
turbed by this accretion process. 



5 Skymaps of the anisotropies 

So far we have focused on analysing particles in cones in the x-, y- and z-direction. We will 
now take a more detailed look on the angular dependence of the /3-profiles. To do so we 
distributed 192 points on a sphere using the HEALPIX framework [32], and defined a cone 
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Figure 11. Skymap of j3 of the major merger with /^initial = 0. 

pointing in the direction of each point. The apex angle was still 45°, so the cone angle is 
larger than the size of one pixel. 

Figure 11 shows a Mollweide projection (which is heavily used in CMB analysis) of /3 
at a radius of 0.35 for the major mergers with /initial = 0. Figure 12 shows the same plot 
for the major merger with an impact parameter along the y-axis, and Figure 13 shows the 
minor merger without an impact parameter. 

In the Mollweide projection the positive x-direction is in the center of the plot, and 
negative x-direction are the points most to the right and left. The positive z-direction is the 
top point, and the negative z-direction is the bottom point. The y-axes are in the two points 
between the positive and negative rc-axis (the positive y-axis is to the right, the negative to 
the left). 

A visible effect is that the presence of an impact parameter breaks the symmetry along 
the y-z plane in Figure 12. It is also clear that the minor merger, mostly affects /3 along the 
infall direction (Figure 13). 

6 /3-7 relations 

We will now compare the merger remnants with the fi-j relation [17], (3 = —0.2 x (7 + 0.8), 
and the attractor [18, 19]. Figure 14 (left panel) shows the spherically averaged f3(^) profiles 
for the major merger remnants. The central panel and the right panel show the same for the 
cones along the x- and y-axis, respectively. 

The /3~7 relations are clearly not obeyed in the two plotted cones, but the spherically 
averaged /3-profiles are in good agreement with the two predictions in the inner parts with 

7 > —2.2. The outer parts with 7 < —2.2 deviate from the relations. It is seen that the 
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Figure 14. f3(j) profiles for the remnant of the major merger between the initial = haloes (with 
and without an impact parameter), and the major merger with two Osipkov-Merritt haloes. The 
remnants are compared between two proposed j3 (7) -relations: f3 = —0.2 x (7 + 0.8) [17] and the 
attractor (the grey lines are taken from [18]). The left panel shows spherically averaged profiles, the 
central panel shows a cone along the collision axis, and right panel shows a cone perpendicular to the 
collision axis. 

spherically averaged /3- and 7-profiles obey the /3-7-relations at the radii, where the /3-profile 
is direction-dependent. 

We see that the spherically averaged ,$(7) profile in the outer parts is strongly dependent 
on a remnants merging history. This finding can explain the large scatter of (3(j) between 
haloes in cosmological simulations [10]. In Figure 2, 3 and 6, we also find that /3(r) depends 
on the detailed merging history, so we conclude that the differences in /3(r) from halo to halo 
in cosmological simulations can be caused by their different merging histories. 

7 Pseudo-phase-space density profiles 

The pseudo-phase-space density profiles (PPSD's), p/o"j? ad and p/a 3 , of the major merger 
remnants with initial = are shown in Figure 15. The profiles are scaled with a factor of 
r 1 j which roughly would give a constant in cosmological simulations [10]. In both remnants 
the normalizations of the PPSD-profiles depend on the cone-direction. The best-fitting value 
for a only has a small variation from cone to cone. 

In several studies dynamical constraints of dark matter haloes have been derived by 
assuming a radial power law behaviour of /o/cj? ad or p/a 3 [33-35]. In the merger remnants 
this assumption is however not correct, due to the different normalizations from cone to cone, 
so such an approximation does not describe the full dynamics of the remnants. 

8 The shape of the haloes 

Is is established that a merger with two spherical haloes results in a triaxial halo with a 
major axis along the collision axis [36]. It has also been shown that the potential and 
velocity contours of a halo are more spherical than the density contours [37, 38] . 

In our merger remnants we measure the minor to major axis ratio, c/a, by comparing 
density profiles in cones through points distributed uniformly of a sphere (with the same 
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Figure 15. The pseudo-phase-space density profiles, /o/cj? ad and p/cr 3 , for the major merger remnant 
with an impact parameter (left panel) and without and impact parameter (right panel). p/c 3 ad is 
multiplied with a factor of 2 for presentation reasons. 



method as in Section 5). To be explicit we will find the maximum and minimum values of 
the radius for a given density, and then calculate the ratio between these two radii. The 
shape will only be determined at r = 1. 

For all the mergers without an impact parameter the halo shapes were aligned with the 
collision axes. For the major mergers with /initial = the axis ratio was 0.59, for the minor 
merger we found a ratio of 0.84, and for the Osipkov-Merritt model we found 0.53. For the 
major merger with an impact parameter we found a ratio of 0.7, and for the minor merger 
with an impact parameter we found 0.86. 

Our simulations are therefore consistent with the discovery in [20], in which it was 
discovered that /3 is aligned with the shape of cosmological haloes, with positive /3-values 
along the major axis, and negative /3-values along the minor axis. 



9 Conclusion 



We have analysed the velocity distributions in different cones centered on merger remnants. 
The velocity anisotropy profiles and the velocity dispersions in the different cones behaved 
differently. This was both the case for major mergers with and without impact parameters, 
and for minor mergers. We also demonstrated that these asymmetries are not washed out by 
a process that mimics smooth accretion. Since mergers are frequent in the real universe, it 
is therefore not surprising that similar asymmetries are present in cosmological haloes [20]. 

In the various merger simulations, the behaviour of (3(r) in the outer parts had a huge 
dependency on the initial conditions. We therefore conclude, that the different behaviour of 
/3(r) from halo to halo in cosmological simulations can be caused by their different merger 
histories. 
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Several studies [39-46] have attempted to derive or characterise distribution functions of 
completely relaxed haloes from first principles. We have demonstrated that merger remnants 
have the merger history encoded in their velocity anisotropy profiles, and they are therefore 
not expected to follow simple distribution functions, where (3 is constant within spherical 
bins or along the isodensity contours. Furthermore it is also clear that the (3-j relation or 
the attractor are not obeyed in all cones in merger remnants. We do, however, find that the 
spherically averaged properties of haloes obey these relations in the inner parts. 

Finally, we note that the merger history likely is important for observational aspects of 
the dark matter haloes. This is exemplified through a correlation between the asymmetric 
nature of the velocity anisotropy of haloes and the surrounding large-scale structure [47], 
along which matter accretion typically occurs [48], and the differences in the line-of-sight 
velocity dispersions of galaxies along the major and minor axes of galaxy clusters [49]. 
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